Observability of Quantum Criticality and a Continuous Supersolid in Atomic Gases 



S. Diehl, 1 ' 2 M. Baranov, 1 ' 2 ' 3 A. J. Daley, 1 ' 2 and P. Zoller 1 ' 2 

1 Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria 
2 Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria 
3 RRC "Kurchatov Institute", Kurchatov Square 1, 123182 Moscow, Russia 

We analyze the Bose-Hubbard model with a three-body hardcore constraint by mapping the 
system to a theory of two coupled bosonic degrees of freedom. We find striking features that could 
be observable in experiments, including a quantum Ising critical point on the transition from atomic 
to dimer superfluidity at unit filling, and a continuous supersolid phase for strongly bound dimers. 



Experiments with atomic quantum degenerate gases 
representing strongly interacting systems have reached a 
level of precision where quantitative tests of elaborate 
many-body theories have become possible [THS]. In the 
interplay between experiment and theory, the challenge 
is now to identify realistic models where quantum fluc- 
tuations lead to qualitatively new features beyond mean 
field in quantum phases and phase transitions. We study 
below the Bose-Hubbard model with a three-body con- 
straint, which arises naturally due to a dynamic suppres- 
sion of three-body loss of atoms occupying a single lattice 
site [6j [7] , and can also be engineered via other methods 
[8] . This constraint stabilizes the system when two-body 
interactions are attractive, allowing for the formation of 
dimers - bound states of two atoms. The phase diagram 
then contains a dimer superfluid (DSF) phase connected 
to a conventional atomic superfluid (ASF). Remarkably, 
this simple but realistic model shows several nongeneric 
features, which are uniquely tied to the three-body con- 
straint and could be observed with cold gases: (i) Emer- 
gence of an Ising quantum critical point (QCP) on the 
ASF-DSF phase transition line as a function of density - 
which generically is preempted by the Coleman- Weinberg 
mechanism [9], where quantum fluctuations drive the 
phase transition first order [T0UT2] , with a finite correla- 
tion length; and (ii) A bicritical point [13] in the strongly 
correlated regime, which is characterized by energetically 
degenerate orders, in our case coexistence of superfluid- 
ity and a charge-density wave, representing a "continuous 
supersolid" with clear experimental signatures. 

Below we describe the constrained model, and then 
present a new analytical formalism for a unified treat- 
ment of onsite constraints in bosonic lattice models, 
based on an exact requantization of the Gutzwiller mean 
field theory. This allows for an analytical treatment of 
the phenomena arising here. At the end we discuss the 
experimental signatures of the latter. 

Constrained model - We consider the Bose-Hubbard 
model on a d-dimensional cubic lattice with a three-body 
onsite hardcore constraint, 

H=-J^ j a\a j -ii^fii + \U y^ni(hj-l), aj~ 3 = 0, (1) 
where denotes summation over nearest neighbors, J 



is the hopping matrix element, \i the chemical potential, 
and U the onsite two-body interaction. The three-body 
constraint stabilizes the attractive bosonic many-body 
system with U < 0, which we focus on here. 

The phase diagram of this model is shown in Fig. [T] 
The dominant phases are an ASF with order parame- 
ters (ai) ^ and (a 2 ) ^ 0, and a DSF with (a 2 ) ^ but 
(di) = 0, formed at sufficiently strong interatomic attrac- 
tion U. In the Gutzwiller mean field approximation [6] 
these two phases are separated by a second order phase 
transition of a special type, the Ising or ASF-DSF tran- 
sition [12] , at which the discrete symmetry of DSF 
is spontaneously broken. One reason to question the 
mean field approach is the presence of two interacting 
soft modes close to the phase transition: the noncritical 
Goldstone mode, related to the (a 2 ) order parameter, 
and the critical Ising mode signaling the onset of atomic 
superfluidity with (a$) ^ 0. This motivates the devel- 
opment of a fully quantum mechanical approach to the 
constrained Hamiltonian 0, with a unified description 
of interaction effects at various scales. 

Formalism: Mapping the constrained model to a cou- 
pled boson theory - Because of the constraint, the op- 
erators ai,a\ are no longer standard bosonic ones and 
the onsite Hilbert space is reduced to only three states 
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FIG. 1: Phase diagram for the Bose-Hubbard model with a 
three-body hard-core constraint, and U < 0. The black curve 
represents the mean field phase border, while red (light gray) 
and blue (dark gray) curves include shifts due to quantum 
fluctuations in d = 2, 3. 
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|a), corresponding to zero (a — 0), single (a = 1), and 
double (a = 2) occupancy. Following Altman and Auer- 
bach [14] , for each lattice site i we introduce three opera- 
tors t\ i creating these states out of an auxiliary vacuum 

state, \a)i = t^Jvac) = (a!) _1//2 (a\^ |vac). By con- 
struction, these operators obey a holonomic constraint, 
^2 a 4 i^oc^i — 1- The Hamiltonian then reads 

where = ^ a ,i- This form is closely related to res- 
onant Feshbach models used, e.g., to describe the BCS- 
BEC crossover in free space [T5] : The "Feshbach term" 
(second line) allows for interconversion of a "dimer" (£2,*) 
into two "atoms" {t\^) on nearby sites, and the detuning 
(first term in the last line), gives the energy difference 
of atoms and dimers. Here the role of the detuning is 
played by the onsite interaction U, in contrast to ~ 1/U 
in the resonant models. 

Using the constraint, we can eliminate one of the op- 
erators t a ^ in Eq. pi, say to,i, as to,i ~^ \fX~i (and 
t\ i — > y/Xl), where Xi = 1 — — hi^. Noting that 
Xf = Xi on the constrained space, we replace \fXi with 
Xi, yielding a polynomial Hamiltonian. The remain- 
ing operators £1,^2 can now be interpreted as standard 
bosonic ones. To demonstrate this, one divides the stan- 
dard bosonic Hilbert space into a physical (Vi) and an 
unphysical (Ui) subspace, Hi = Vi®Ui, where the phys- 
ical one consists of states obeying the 3-body constraint. 
We see that H does not couple the physical V = \\Vi 
and unphysical U = Y[U{ subspaces. 

The distinction between the contributions from phys- 
ical and unphysical spaces can most conveniently be 
achieved by using the quantum effective action r[ti,t2] 
[T6] , which is a functional on classical fields and provides 
all one-particle irreducible correlation functions as coeffi- 
cients of an expansion in powers of t\ and t^. Because T is 
formulated in terms of physical quantities, we can restrict 
its general form to a polynomial obeying the three-body 
constraint. The form of the effective action thus is re- 
stricted by the constraint, in addition to the symmetries 
of the microscopic theory. 

To apply the above construction to a many-body sys- 
tem, one has to choose the proper ground state and the 
corresponding operators. Following Ref. [17], we intro- 
duce new operators b a ^ = (Ri) a ptp,i (<^,/3 = 0,1,2), 
with a unitary transformation R. The parameters of R 
are such that 60, % creates the mean field vacuum, and 
bi y i and 62, i correspond to fluctuations on top of this 
state, with vanishing expectation values (see [15]). The 
DSF ground state, for example, corresponds to: bo,i = 



cos(#/2)£ ,i + sin(0/2)exp(-z0)£ 2 ,i, h,i = cos(#/2)£ 2 ,i - 
sin(0/2) exp(i0)t o ,i, and bi yi = where <p is an ar- 
bitrary phase and the angle G [0,7r] is such that 
2sin 2 (#/2) = n, the density of atoms (on the mean field 
level for simplicity). The operators are subject to 
the same constraint, b^ a — 1, and we can elim- 
inate b^l as described above. This results in a polyno- 
mial Hamiltonian for the remaining operators b\ ^ and 
&2,i, where the operator independent part precisely re- 
produces the Gutzwiller mean field energy. 

Application of the formalism - The limit n — > gives 
a stringent check of our formalism, where we recover 
the nonperturbative Schrodinger equation for the dimer 
bound state formation (see [18] for details). We now ap- 
ply our method in the many-body context: 

(i) Ising quantum critical point. The polynomial 
Hamiltonian describes atomic (bi,$) and dimer (&2,i) fluc- 
tuations around the spatially uniform DSF state. Af- 
ter taking the long wavelength continuum limit, one can 
easily see that there are two soft modes in the vicinity of 
the ASF-DSF transition: the noncritical Goldstone mode 
7r ~ 1111(62) corresponding to the U(\)/Z2 gauge sym- 
metry broken by the presence of the dimer condensate, 
and the critical atomic Ising mode cp ~ Re(&i) signaling 
the appearance of an atomic condensate. The other two 
modes remain massive and do not affect the physics of 
the phase transition. Integrating out the latter we obtain 
an effective low energy action for the soft modes [T5] : 

S e ^7r} = f {^(-Z^-^A + m^ + A^ 4 

J X 

+ |tt(-Z^ - £ 2 A)tt + i^ 2 <9 r 7r}. (3) 

It describes phonons 7r in the dimer superfluid coupled 
to a real scalar Ising field tp, in turn represented by an 
action of the Ginzburg-Landau type with the "mass" 
parameter crossing zero at the ASF-DSF transi- 
tion. The coupling n comes from the cubic coupling 
— ^/2J cos(6)b\ jbi^bij + h.c. which originated from the 
"Feshbach term" in the original Hamiltonian Eq. ([2|, 
such that k ~ cos(#) « 1 — n. This cubic coupling of 
Goldstone to Ising mode with linear time derivative has 
the same degree of relevance as the Ising coupling A, lead- 
ing to the Coleman- Weinberg (CW) phenomenon [TO] : 
The renormalized value of the Ising coupling A reaches 
zero at some finite scale £ at which m\ is still positive. 
As a result, terms with higher powers of (p, which are 
generated by fluctuations, become important. These self- 
interaction terms provide a new minimum with (cp) 7^ 0, 
which is reached via a first order phase transition with fi- 
nite correlation length £. Therefore, the ASF-DSF phase 
transition in our model is actually first order, contrary 
to the predictions of the mean field approach. 

The coupling via a temporal derivative and, therefore, 
the CW phenomenon, is rather generic in nonrelativistic 
systems, in which an Ising field emerges as an effective 
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low energy degree of freedom [T0HT2] . In our case, how- 
ever, the coupling k vanishes for n —> 1. The existence of 
such a decoupling point can be proven assuming a contin- 
uous, monotonic behavior of a particular compressibility, 
the \i -derivative of the dimer mass term K = —dm 2 d /dn\ n 
within the DSF phase: it then must have a unique zero 
crossing, because it is > (<)0 for n = 0(2) [18]. This 
argument is tied to the existence of a maximum filling, 
and thus to the three-body constraint. Using the Ward 
identities resulting from a temporally local gauge invari- 
ance b a — » exp(iaA(£))6 a , a = 1,2, and \i — >• \i + iA(t) 
[T9] , we see that k oc K. As a result, the first order Z 2 
transition terminates into a true d + 1 dimensional Ising 
quantum critical point in the vicinity of unit filling. 

Fluctuations also shift the mean field phase boundary, 
cf. Fig. [T] The shift is only pronounced for n <C 1, the 
reason being that dimer formation and atom criticality 
approach each other for n — » [T8] . 

(ii) Continuous supersolid. Another peculiar conse- 
quence of the three-body constraint occurs for dimers 
in the strong coupling limit, where single particle exci- 
tations are strongly gapped (~ \U\/2) and can be inte- 
grated out perturbatively. Taking the dominant nearest 
neighbor hopping t and interaction v into account, the 
effective lattice theory for hardcore bosons (dimers or 
di-holes) can be conveniently rewritten as an antiferro- 
magnetic Heisenberg spin Hamiltonian (see, e.g., [21]), 

Haf = 2t ^2 I s 

Syj+\s Zii s Zij ] (4) 

(id) 

restricted to a subspace with a fixed projection of the 
total spin on the z-axis, S z = J2i s z,i — L{na — 1/2), 
where L is the total number of the lattice sites, and 
nd = n/2 the dimer filling. The anisotropy parameter 
A is the ratio of the interaction and hopping, A = v/2t. 
In the leading second order perturbation theory, we find 
t = v/2 = 2J 2 /\U\, such that A = 1 and the Hamiltonian 
Q is SO (3) -invariant, corresponding to a symmetry en- 
hancement compared to the conventional SO (2) ~ U(l) 
phase symmetry for bosons. It parallels a similar ef- 
fect for attractive lattice fermions [22], and is a pecu- 
liar feature of the three-body hardcore constraint - if 
virtual triple and higher occupancies are allowed, one 
finds A = 4 [23 . The symmetry enhancement is oper- 
ative for exactly half filling of dimers, rid = 1/2, where 
S z = 0, while for other dimer fillings S z ^ and the 
symmetry is reduced to U(l). In the first case, however, 
the ground state of the Hamiltonian Q is the antiferro- 
magnetic state parametrized by the direction of the Neel 
order parameter on the three-dimensional Bloch sphere. 
Generically, the Neel order parameter has components 
both in the xy -plane and along the z -axis. This means 
that the ground state of bosons has both DSF and charge- 
density wave (checkerboard- like, CDW) orders, i.e. is 
a supersolid. The specific feature, however, is that the 
SO (3) symmetry admits a continuous change in the ratio 



between the two order parameters without changing the 
energy, and this state may thus be called a continuous 
supersolid. A particular choice of the order parameters 
depends on the way the system is prepared and on the 
boundary conditions. This is in contrast to other occur- 
rences of supersolidity in bosonic systems [24] . 

The spontaneously broken SO (3) symmetry in the con- 
tinuous supersolid provides us with two massless Gold- 
stone modes. A spin wave analysis yields their dispersion 

^(q)=^[? q (l + A-A? q )] 1/2 

with e q = l/d^2 x (l — cosqe A ), A < 1. For A = 1 the 
second Goldstone mode emerges at the edge of the Bril- 
louin zone, adding to the one at zero momentum. In 
the next (fourth) order of perturbation theory, we find 
[TBI A = 1 - S(z - 1)(J/ \U\) 2 < 1, and the DSF ground 
state is slightly favored over CDW order due to weak 
explicit breaking of SO '(3). Still, the proximity to the 
continuous supersolid manifests itself in a weakly gapped 
(A ~ tz(l — A)) collective mode at the edge of the Bril- 
louin zone, which may be probed experimentally, see be- 
low. 

Experimental signatures - We now discuss in detail the 
experimental signatures of the critical behavior and con- 
tinuous supersolid. Above we showed that the ASF-DSF 
phase transition (see Fig. 1) terminates in an Ising quan- 
tum critical point at n = 1. The phase transition is thus 
second order at n = 1, with a diverging correlation length 
at the transition point. Away from n = 1 the transition 
is first order, with a finite correlation length estimated as 
£/a ~ k~ 6 ~ |1 — n\~ 6 (a the lattice spacing) using the 
renormalization group flow of [10 . As typical of the ra- 
diatively induced first order transitions, the near-critical 
domain is large, with a correlation length exceeding the 
typical extent of optical lattices of 50 to 100 sites in a 
region 1/2 < n < 3/2. This behavior can be measured 
directly in experiments probing the correlation length, 
as done in Ref. [1] . Alternatively, critical opalescence via 
damping of collective oscillations [20 could be probed. 

The continuous supersolid appearing at strong attrac- 
tion and unit filling can be detected by measurement of 
the coexisting spatial order and dimer superfluid correla- 
tions. The spatial structure could be detected via noise 
correlation measurements [25], and the dimer superfluid 
correlation functions in the momentum distribution of 
dimers, which could be imaged, e.g., by associating atoms 
to molecules, and measuring their momentum distribu- 
tion. Strong evidence would also be obtained by probing 
collective modes that appear at the edge of the Brillouin 
zone. It is also possible to stabilize the CDW by ramping 
a weak (~ A <C tz) superlattice, acting as a staggered 
external field that rotates the Neel order parameter from 
the xy -plane to the z -axis. 

We further elaborate on experimental signatures using 
exact numerical time-dependent density matrix renor- 
malization group (t-DMRG) techniques [26] in ID in the 
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FIG. 2: Ground state computations on 60 lattice sites with 
JJ I J — —20 using t-DMRG. (a) Correlation functions charac- 
terizing the CDW and DSF orders for open boundary condi- 
tions at unit filling on a log-log scale as a function of distance 
x. (b) Fitted algebraic decay exponents Kdsf and i^cDW for 
varying n (errorbars show estimates of the fitting error), (c) 
Density n(x) for a system with a harmonic trapping potential 
V(x)/J = {x- 30.5) 2 /900, N = 60 particles, (d) Shaded plot 
of the DSF correlation function {bib^byby) with interpolated 
shading, indicating substantial DSF order where n ~ 1 . 



ground state for realistic experimental size scales and pa- 
rameters: Fig. 2a shows the density-density correlation 
function (n^n^) — (ni)(rii+ x ) characterizing CDW or- 
der, and the DSF correlator (b\blbi+ x bi+ x ) as a function 
of distance x. At unit filling, both decay algebraically 
and are essentially equal, indicating coexistence both or- 
ders. In Fig. 2b we show the result of fitting an algebraic 
decay x XcDW DSF to the correlation functions. For n = 1 
these are equal, but away from unit filling, the DSF cor- 
relations decay more slowly, so that DSF order dominates 
CDW 1 . In experiments in a harmonic trap, where the 
filling factor varies across the system, this gives rise to 
further signatures. In Fig. 2c, we plot the density in a 
trap, showing that a region exists near unit filling where 
oscillations in the density are present, characteristic of 
the appearance of CDW order. Fig. 2d shows that in the 
same region the DSF correlations are significant, whereas 
in the center of the trap, a constraint-induced insulating 
phase with n — 2 appears. For more details see [T8] . 

Conclusion - We have demonstrated that an atomic 
Bose gas in an optical lattice with three-body onsite con- 
straint provides a realization of such fundamental phys- 
ical concepts as the Coleman- Weinberg phenomenon of 
radiative mass generation and Ising quantum criticality. 
In addition, the ground state at unit filling in the strongly 



correlated limit is an example of a continuous supersolid 
- a supersolid with a tunable ratio between the superfluid 
and the charge-density wave order parameters. 
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